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1.  INTRODUCTION 


1.1  Overview 

Science  Research  Laboratory  has  developed  a  family  of  all- solid-state  pulsed  power  gen¬ 
erators  featuring  levels  of  efficiency,  reliability  and  power  density  far  above  what  had  been 
previously  obtainable.  A  common  feature  of  these  systems  is  SCR-commutated  nonlinear  mag¬ 
netic  pulse  compressors.  Using  SCRs  rather  than  the  more  traditional  thyratrons  decreases 
system  weight  and  volume  while  providing  a  significant  increase  in  reliability.  Pulsers  using  this 
technology  have  demonstrated  lifetimes  of  10*°  shots. 

The  primary  objective  of  this  SBIR  is  to  find  ways  of  combining  these  pulsed  power  systems 
with  new  forms  of  CO2  laser  amplifiers  as  the  first  step  in  the  development  of  a  new  family  of 
COi  laser  radar  systems  suitable  for  space-based  applications.  The  CO2  lasers  investigated  in 
this  SBIR  feature  conduction-cooling  of  the  laser  mixture  rather  than  the  more  typical  flow  loops 
and  heat  exchangers.  Conduction  cooling  has  significant  size,  weight  and  system  complexity 
advantages  which  complement  those  of  the  all-solid-state  pulsed  power. 


1.2  Program  Plan 

In  this  SBIR  we  built  two  forms  of  conduction-cooled  CO2  lasers  matched  to  separate 
forms  of  SRL’s  solid-state  pulsers.  The  first  system  uses  a  spark  gap  preionized,  transverse 
discharge,  rectangular  slab  geometiy  laser  amplifier  matched  to  an  existing  solid-state  pulser 
(COLD-I).  The  laser  medium  in  this  slab  system  was  interferometrically  probed  with  a  pulsed 
laser  to  test  the  effectiveness  of  its  conduction  cooling.  The  second  system  uses  a  longitudinal 
discharge,  cylindrical  geometry  laser  amplifier  powered  by  a  new  pulsed  power  system  (SSLAM- 
Vm)  capable  of  much  higher  output  power.  This  new  pulsed  power  system  features  separate 
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modules  which  individually  produce  12  kV,  5  J  pulses  but  which  can  be  stacked  to  provide  the 
72  kV,  30  J  pulses  required  for  the  longitudinal  discharge.  The  ouq)ut  power  of  this  second 
system  was  optimized  by  varying  laser  gas  mixtures,  pressure  and  pulse  repetition  rate. 

U  Summary  of  Results 

The  interferometric  measurements  on  the  transverse  discharge,  slab  geometry  lasers  showed 
that  the  thermal  energy  deposited  by  a  single  pulse  dissipated  completely  within  100  ms.  This 
means  that,  with  actively  cooled  walls,  the  system  could  be  run  continuously  at  10  Hz  without 
the  heat  from  one  pulse  affecting  the  next.  At  50  ms  the  interferometric  path  error  was  only 
a  small  fraction  of  a  visible  wave.  The  test  fringes  showed  the  characteristic  cylindrical  error 
expected  of  a  slab  which  is  heated  uniformly  and  cooled  at  its  walls.  This  indicates  that  the 
system  could  have  been  run  at  20  Hz  with  only  a  small  effect  on  laser  beam  quality.  The  existing 
pulsed  power  supply  repetition  rate,  however,  was  limited  to  10  Hz. 

The  longitudinal  discharge,  cylindrical  geometry  laser  was  run  successfully  at  repetition 
rates  of  up  to  100  Hz.  Its  average  output  power  typically  peaked  at  50  Hz,  depending  on  the 
gas  mixture  and  pressure.  The  50  Hz  practical  limit  was  expected  as  calculations  had  indicated 
that  the  core  of  the  discharge  would  be  heated  excessively  above  that  rate.  The  efficiency 
of  conversion  from  energy  deposited  in  the  discharge  to  laser  output  energy  also  varied  with 
gas  mixture  and  pressure  but  could  be  brought  up  to  about  20%.  The  maximum  pulse  energy 
obtained,  however,  was  only  350  mJ  for  a  maximum  average  power  of  18  W.  The  essential 
problem  is  that  the  optimum  load  for  the  SSLAM-VHI  pulser  was  measured  to  be  about  4  12 
whereas  the  typical  discharge  impedance  ran  in  the  200  to  400  f2  range.  A  likely  solution  to 
this  problem  has  been  suggested  but  not  yet  implemented. 
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2.  TRANS^TRSE  DISCHARGE  SYSTEM  TESTS 


2.1  Laser  Design  Summary 

While  the  new  SSLAM-VIII  pulsed  power  system  was  being  built  for  this  project,  tests 
of  conduction  cooling  were  performed  on  a  different  CO2  laser  prototype.  These  experiments 
utilized  a  new  version  of  a  conduction-cooled,  transverse  discharge,  slab  geometry  CO2  laser  sys¬ 
tem.  An  earlier  version  of  this  system  had  been  produced  under  a  previous  contract^*)  but  proper 
tests  of  its  medium  quality  had  not  been  completed.  This  report  will  describe  interferometric 
tests  of  the  laser  medium  quality  for  the  new  version  of  the  slab  geometry,  conduction-cooled 
laser. 

Both  lasers  utilize  a  cylindrical  housing  to  contain  both  the  discharge  electrodes  and  the 
low  pressure  (30  Torr)  laser  gas  mixture  required  for  conduction  cooling.  Two  glass  plate 
walls  separate  the  discharge  electrodes  and  serve  to  constrain  the  discharge  within  the  enclosed, 
rectangular  slab  volume.  In  both  systems  the  dimensions  of  the  active  volume  were  13.5  cm 
in  the  direction  of  the  discharge,  2.4  cm  between  the  two  glass  plates  and  100  cm  along  the 
optical  cavity  axis  (a  total  of  3.2  liters).  The  distinctive  slab  geometry  of  this  laser  results  from 
the  requirement  that  heat  produced  in  the  discharge  be  carried  by  conduction  to  the  laser  walls. 
Figure  1  shows  a  photograph  of  the  discharge  produced  by  the  new  system  along  with  a  number 
of  the  discharge  parameters.  If  the  walls  are  too  far  apart  the  heat  cannot  be  dissipated  before  the 
next  pulse.  If  the  walls  are  too  close  the  total  gas  volume  is  insufficient  to  provide  the  required 
energy  per  pulse.  A  detailed  discussion  of  the  constraints  for  a  conduction-cooled  system  are 
given  in  the  proposal  for  this  Phase  11  SBIR.^^) 

The  new  version  of  the  system  employed  several  improvements  over  the  original  design. 
The  first  version  of  the  laser  had  utilized  a  large  number  of  individual  cables  feeding  through 
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CONDUCTION  COOLED  CO2  LASER  DISCHARGE 


\'oluin<‘ 

\ 

(  t Kxlc  Spiu'ing 

Cm  H  ilt 

Discliait',!'  Pulse  L(Mip;th 
Lasci  Pulse  Length 

(  )VIM  ,lll  fvffit  il'IM  V 
Total  (Jas  Ph'ssuh’ 


2.4  X  13.5  X  100  env' 

5  kV 

13.5  (  111 

SOO  Ainp.s 

3  fis 

40  /IS 

30  Tori 


Figure  1 
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the  cylinder  wall  to  power  the  spark  gap  preionizers  spaced  along  one  of  the  electrodes.  Those 
cables  tended  to  leak  which  compromised  the  laser  gas  purity.  The  new  version  of  the  system 
utilized  only  two  cable  feedthroughs  and  each  spark  pin  is  individually  ballasted. 

2.2  Interferometric  Probe  Design 

The  physics  of  the  conduction  cooling  process  is  sufficiently  straightforward  that  its  lim¬ 
itations  were  well  understood.  The  interferometric  tests  on  the  laser  medium  were  primarily 
designed  to  address  the  question  of  whether  or  not  there  would  be  high  order  distortions  due 
to  residual  sound  waves  or  shock  waves  in  the  medium  left  over  from  previous  shots.  An 
interferometric  probe  was  built  to  test  for  these  distortions. 

The  interferometric  diagnostic  consists  of  a  Mach-Zehnder  interferometer  and  a  spatially 
filtered,  pulsed  dye  laser  probe  beam.  The  dye  laser  pulse  length  is  only  3  ns  so  the  interferometric 
probe  is  a  “snapshot”  of  the  state  of  the  laser  medium  at  the  time  of  the  probe.  A  digital  delay 
generator  was  used  to  sweep  the  time  of  the  probe  laser  relative  to  the  CO2  laser  discharge 
on  successive  shots.  A  CCD  camera  and  digital  framegrabber  was  used  to  record  the  fringes 
produced.  The  wavefront  which  bests  fits  these  fringes  was  determined  with  an  algorithm 
developed  at  SRL.  This  algorithm,  presented  in  the  Appendix,  calculates  a  wavefront  defined 
over  an  NxN  grid.  A  minimum  surface  area  constraint  on  the  fit  wavefront  defines  those  grid 
values  which  are  not  directly  determined  from  the  fringes. 

23  Interferometric  Probe  Data 

Figure  2  shows  a  partial  set  of  the  fringe  data  taken.  Figures  3  to  7  show  the  wavefront 
obtained  from  these  firinges  using  an  iterative,  grid  fit  algorithm  developed  by  SRL  (listed  in  the 
Appendix). 

The  top  graph  in  Fig.  3  shows  the  fringes  taken  just  before  a  discharge  pulse  with  a  wall 
temperature  roughly  equal  to  the  initial  gas  temperature.  It  takes  a  few  microseconds  for  the  gas 
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Figure  2:  Interferometric  probes  of  the  transverse  discharge,  conduction 
cooled,  CO2  laser  medium  taken  at  the  spark  preionizer  end 
of  the  discharce. 
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heated  by  the  discharge  to  expand.  The  optical  path  difference  (OPD)  thus  does  not  reach  its 
maximum  value  for  about  SO  fis.  In  the  1  ms  time  range  the  OPD  is  complicated  by  what  appears 
to  be  cooling  in  the  center  of  the  discharge  channel.  This  effect  is  likely  caused  by  energy  stored 
in  vibrational  modes  of  the  gases  involved  and,  although  interesting,  does  not  affect  the  laser  in 
any  practical  way.  The  OPD  at  20  ms  shown  in  Fig.  6  is  the  error  which  the  laser  pulse  would 
see  on  the  next  pulse  at  a  SO  Hz  repetition  rate.  The  error  of  successive  pulse  add  so  that  the 
error  on  the  third  pulse  at  SO  Hz  would  be  about  the  error  at  20  ms  plus  the  error  at  SO  ms.  The 
total  error  at  SO  Hz  would  thus  still  be  below  1/10  of  a  wave  at  the  CO2  wavelength.  The  slight 
rise  in  temperature  in  the  center  of  the  cavity  will  cause  the  discharge  to  strike  preferentially  at 
that  point.  This  feedback  effect  was  noted  in  the  longitudinal  laser  system  to  be  described  in  the 
next  section. 

The  transverse  discharge  slab  laser  could  only  be  fired  at  a  maximum  rate  of  10  Hz  due 
to  the  choice  of  capacitors  used  in  COLD-I.  The  OPD  error  at  95  ms  shown  in  Fig.  7  was 
essentially  at  or  below  the  measurement  limit  of  the  interferometer  used.  This  means  that  the 
transverse  discharge  system  always  appeared  to  be  running  in  the  single  pulse  mode  as  far  as 
the  laser  medium  quality  was  concerned. 

The  experimental  plan  was  to  continue  the  measurements  on  the  longitudinal  discharge, 
conduction  cooled  laser  system  which  could  be  run  at  100  Hz  and  above.  In  practice,  however, 
the  efficiency  of  the  energy  deposition  in  the  longitudinal  discharge  was  too  low  to  provide  a 
good  test  of  the  conduction  cooling  in  that  geometry. 
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3.  LONGITUDINAL  DISCHARGE  SYSTEM  TESTS 


3.1  SSLAM-Vm  Pulsed  Power  System 

The  decision  to  develop  a  longitudinal  discharge  laser  system  was  prompted  by  the  sig¬ 
nificant  advantages  of  the  SSLAM  series  of  pulsed  power  generators  over  more  conventional 
systems.  These  pulsed  power  supplies  are  nearly  indestructible,  require  an  input  voltage  of  only 
600  V  and  can  operate  continuously  at  5  kHz.  In  the  SSLAM  series  a  high  voltage  is  induced  on 
an  isolated  conducting  rod  using  techniques  borrowed  from  induction  accelerators.  The  resulting 
power  supply  can  be  built  in  a  modular  fashion  with  successive  modules  simply  adding  another 
step  in  voltage.  A  longitudinal  discharge  requires  these  high  voltages  and  is  thus  an  obvious 
candidate  for  the  application  of  this  system. 

Figure  8  shows  a  single  SSLAM-Vin  module  which  can  produce  5  J  pulses  at  12  kV  at 
repetition  rates  of  up  to  5  kHz.  The  module  measures  7  inches  in  height  by  17  inches  in  width 
and  28  inches  in  depth  and  weighs  39  kg.  Figure  9  shows  a  close  up  of  the  two  induction  rods 
in  each  module  and  how  they  are  connected  in  series  to  produce  a  total  of  12  kV.  Figure  10 
shows  the  rack  panel  containing  six  12  kV  modules  which  produces  a  total  pulse  energy  of  30 
J  at  72  kV. 

32  Laser  Design  Summary 

To  match  the  high  voltages  which  could  be  obtained  from  SSLAM-Vm,  a  longitudinal 
discharge  laser  was  designed.  The  discharge  cell  consists  of  a  100  cm  long  thin  wall  glass  tube 
with  a  4.13  cm  inside  diameter.  The  voltages  available  from  SSLAM- Vm  are  sufficient  to  break 
down  the  gas  at  pressures  of  up  to  about  30  Torr,  which  are  also  the  pressures  appropriate  for 
conduction  cooling. 
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The  ends  of  the  tube  are  fitted  into  cylindrical  electrodes  which  also  serve  to  retain  ZnSe 
windows  to  seal  off  the  system.  The  SSLAM-Vin  power  supply  modules  were  connected  to 
supply  a  potential  of  +35  kV  to  the  electrode  at  one  end  of  the  discharge  and  -35  kV  to  the 
electrode  at  the  other  end.  A  grounded  ring  near  the  electrodes  at  the  ends  of  the  cavity,  but 
outside  of  the  glass  tube,  served  to  break  down  the  gas  near  the  electrodes  at  the  beginning  of 
the  pulse  and  thus  initiate  the  full  discharge. 

The  laser  gas  mixture  was  fed  into  the  laser  cavity  through  the  electrodes  at  one  end  of 
the  cavity  and  extracted  through  a  bleed  valve  at  the  other  end.  Generally  the  laser  mixture  was 
not  circulated  in  these  tests.  The  system  was  also  only  run  in  bursts  of  up  to  several  thousand 
pulses  to  avoid  significant  heating  the  walls  of  the  discharge  tube. 

3J  Conduction  Cooling  in  a  Cylindrical  Geometry 

Since  the  original  design^^)  had  considered  a  slab  geometry  it  is  worthwhile  to  recalculate 
the  effects  of  conduction  cooling  for  this  cylindrical  geometry.  The  general  steady-state  equation 
for  the  temperature  T  in  a  gas  of  conductivity  k  heated  at  a  rate  Q  (power  per  unit  volume)  is 

V.(«VT)  =  -Q  (1) 

In  general,  both  Q  and  «  (which  is  temperature  dependent)  can  be  functions  of  position  in  the 
gas. 

In  a  cylindrically  symmetric  geometry,  Eq.  (1)  reduces  to 


where  r  is  the  radial  position  (0  <  r  <  R). 

Eq.  (2)  can  be  solved  numerically  for  a  given  heating-rate  distribution  and  temperature- 
dependent  conductivity.  A  simplification  occurs  if  both  Q  and  k  are  independent  of  r;  for  this 
case  an  analytical  solution  may  be  obtained  immediately.  One  has 

^^_Qr 

dr  k2 
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Let  Ttt,  =  T(R),  the  wall  temperature.  Then 


T(r)-T«,  =  £(il2-r2)  (3) 

and  thus  the  maximum  temperature  obtained,  Tmax  =  T'(O)  is 

I'maz  =Tv,+  (^) 

where  D  =  2R  is  the  inner  diameter  of  the  tube. 


It  is  interesting  to  compare  this  result  with  the  ma3cimum  temperature  in  a  cell  having  a 
rectangular  cross  section  of  dimensions  d  x  cf  in  the  limit  d  «  ct.  In  that  limit  the  conduction 
to  the  walls  separated  by  d  dominates,  and  the  conduction  to  the  walls  separated  by  cf  may  be 
neglected.  For  the  same  heating  rate  Q  the  result,  easily  obtained  from  Eq.  (1)  is  then 

Tmax  =  Tu,  +  (5) 

Thus  as  far  as  maximum  temperature  increase  is  concerned,  a  tube  of  circular  cross  section  and 
diameter  D  is  equivalent  to  a  tube  of  rectangular  cross  section  and  small-wall  separation  D/y/2. 
This  result  is  reasonable:  the  cylindrical  tube  has  more  nearby  surface  area  through  which  the 
heat  can  flow  than  a  rectangular-cross-sectional  tube  with  the  same  wall  separation. 


Since  the  volume  cf  the  heated  cylinder  is 


V  = 


where  L  is  the  length  of  the  cylinder,  the  maximum  temperature  rise  of  the  cylinder  may  be 
expressed  in  terms  of  the  total  heating  rate 


as  follows; 


Qt  =  VQ  =  irR'^LQ 
Qt 


I'max  — 


(6) 


(7) 


In  addition,  the  total  heating  rate  may  be  related  to  the  laser  power  P  and  efiiciency  e  (defined 
as  the  ratio  of  the  optical  power  extracted  to  the  electrical  power  pumping  the  laser  gas)  by 

(1  -  e)P 


Qt  = 


e 
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(8) 
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so  that 


(9) 


r  r  -(1-^)^ 

^ireKL 

For  a  3:2:1  Hc:N2:C02  mixture  typical  of  a  CO2  laser,  the  conductivity  is  ~  1  x  10“^ 
W/cm  °K.  Thus  for  a  cylindrical  tube  of  length  L,  laser  power  P  and  efficiency  e  the  steady-state 
temperamre  rise  at  the  center  is 

Tmax  -  ^  80^^^  (10) 

where  P  is  measured  in  watts  and  L  is  measured  in  cm. 

Consider  for  example  a  CO2  laser  of  length  L  =  100  cm,  average  power  P  =  20  W  and 
efficiency  e  =  10%.  According  to  Eq.  (10),  the  maximum  time-averaged  temperature  rise  is 

AT  ~  180“K 

Temperatures  on  this  order,  of  course,  will  drop  the  gas  density  on  the  longitudinal  axis  to 
about  0.6  of  the  gas  density  at  the  wall.  The  discharge  will  then  preferentially  pump  the  laser 
mixture  on  the  axis.  This  effect  was  observed  in  the  experiments  described  below. 


3.4  Optimization  of  the  Laser  Output  Power 

Figures  11  to  14  show  the  voltage  and  current  through  the  longitudinal  discharge  at  pulse 
repetition  rates  of  from  20  Hz  to  100  Hz.  In  all  of  these  particular  cases  the  laser  gas  was  1:1:8 
mixture  of  C02:N2:He  and  only  the  data  for  the  30th  pulse  in  a  burst  of  30  is  shown. 

The  voltages  were  measured  with  an  optically  isolated  probe  consisting  of  an  LED  trans¬ 
mitter  attached  between  the  two  discharge  electrodes  in  series  with  a  2  Mfl  current  limiting 
resistor.  The  light  ouqrut  from  the  LED  is  proportional  to  the  current  through  the  LED  and  thus 
to  the  voltage  across  the  discharge.  The  LED  was  coupled  through  an  insulating  optical  fiber  to 
a  fast  photodiode  receiver.  The  entire  probe  was  calibrated  with  fast,  high  voltage  pulses  from  a 
known  source.  The  current  was  measured  with  a  fast  Pearson  current  transformer  which  is  also 
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Discharge  Power  (  MW  )  Voltage  (  kV ) 


Longitudinal  Laser  Discharge,  20  Hz 


150 


100 


50 
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Figure  11:  Voltage  and  current  waveforms  for  a  typical  laser  gas  mixture 
and  pressure  at  a  repetition  rate  of  20  Hz. 
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Average  Input  Power  (  W  ) 


Discharge  Power  (  MW  )  Voltage  (  kV  ) 


Longitudinal  Last^r  Discharge,  30  Hz 


0  12  3  4 


Time  (^is) 

Figure  12:  Voltage  and  current  waveforms  for  a  typical  laser  gas  mixture 
and  pressure  at  a  repetition  rate  of  30  Hz. 


23 


SCIENCE  RESEARCH  LABORATORY 


Average  Input  Power  (  W  ) 


Discharge  Power  (  MW  )  Voltage  (  kV  ) 


Longitudinal  Laser  Discharge,  50  Hz 


0  12  3  4 


Time  ( ) 

Figure  13:  Voltage  ’ind  current  wavefoms  for  a  typical  laser  gas  mixture 
and  pressure  at  a  repetition  rate  of  50  Hz. 
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Average  Input  Power  ( W  ) 
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Longitudinal  Laser  Discharge,  100  Hz 


350 

300 

250 

200 

150 

100 

50 

0 

8 

7 

6  -- 


5  >. 

O) 

0) 

4  c 
^  UJ 

a> 

3  1“ 

JZ 

u 

p  S2 

b 


1 

0 


< 


c 


3 

o 


800 

700 

600 

500 

400 

300 

200 

100 

0 


Figure  14:  Voltage  and  current  waveforms  for  a  typical  laser  gas  mixture 
and  pressure  at  a  repetition  rate  of  100  Hz. 
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Average  Input  Power  (  W  ) 


isolated  from  the  discharge  circuit.  Both  the  current  and  voltage  waveforms  were  digitized  at 
100  MSamples/s  with  a  Lecroy  oscilloscope  and  stored  for  later  processing.  The  input  power 
and  input  energies  shown  on  the  lower  graphs  are  calculated  from  these  voltage  and  current 
waveforms.  A  separate  digital  scope  stored  the  energy  of  the  laser  pulse  as  measured  by  an 
integrating  pyroelectric  joulemeter. 

A  total  of  31  successive  pulse  waveforms  and  energies  were  stored  for  each  run,  starting 
with  the  first  trigger  to  the  SSLAM-Vin  pulser.  The  first  trigger  pulse  does  not  result  in  a 
discharge  since  the  SSLAM-Vin  system  uses  a  resonant  charging  system.  The  first  pulse  is  thus 
only  at  half  voltage  and  the  gas  does  not  break  down.  The  wall  temperatures  before  each  run 
were  nominally  at  room  temperature  and  the  walls  were  sufficiently  thick  that  they  did  not  heat 
significantly  during  the  run. 

At  low  repetition  rates  the  gas  cools  sufficiently  for  each  pulse  to  be  independent,  except 
for  a  slow  energy  decay  due  to  the  dissociation  of  some  of  the  CO2.  Figure  15  shows  the 
energy  for  each  pulse  in  the  30  pulse  series  at  the  repetition  rates  illustrated  above.  At  higher 
repetition  rates  the  behavior  of  the  first  few  pulses  is  generally  erratic,  exhibiting  higher  energies 
but  sometimes  skipping  pulses. 

In  all  of  these  cases  the  voltage  typically  rises  to  20  to  30  kV  before  the  discharge  is 
initiated.  The  current  spikes  for  about  150  ns  and  then  falls  back  to  a  more  or  less  steady 
state  value  for  about  1  fis.  About  half  of  the  energy  transferred  to  the  discharge  occurs  during 
the  initial  voltage  spike  when  the  discharge  impedance  is  below  100  il.  The  impedance  rises 
following  the  initial  current  spike  to  about  1000  Q  and  then  drops  again  to  around  200  Q  for 
most  of  the  rest  of  the  discharge.  Tests  on  the  SSLAM-Vin  pulser  show  that  it  is  best  matched 
to  much  lower  impedances  of  about  4  Q.  This  is  the  reason  that  the  total  energies  transferred 
to  the  discharge  are  typically  only  about  5  to  6  J  rather  than  the  30  J  that  the  pulser  can  deliver 
to  a  matched  load.  It  is  likely  that  the  overall  efficiency  of  the  system  could  be  improved  by 
effectively  pulse  charging  a  local  energy  storage  capacitor  attached  to  the  discharge  tube.  The 
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Energy  (J  ) 


Energy  vs.  Pulse  Position 


0  5  10  15  20  25  30 


Shot  Number  In  Burst 

Figure  IS:  Eneigies  for  each  laser  pulse  resulting  from  31  successive 
triggers  to  the  SSLAM-VHl  pulsed  power  system. 
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charging  would  occur  through  a  string  of  high  voltage  diodes  which  would  block  the  return 
current  when  the  SSLAM-Vin  voltage  pulse  drops  back  to  zero.  We  were  not  able  to  implement 
this  iix  before  the  end  of  the  contract,  however. 

The  voltage  and  current  traces  are  essentially  independent  of  the  repetition  rate  of  the 
discharge.  This  is  not  the  case  for  the  laser  energy.  Pulse  energies  peak  at  the  30  Hz  repetition 
rate  as  shown  in  Fig.  IS.  At  repetition  rates  above  20  Hz  the  discharge  appears  to  stick  close 
to  the  walls  of  the  glass  tube  for  the  first  few  shots  and  then  concentrate  near  the  center  of  the 
tube  after  the  gas  has  heated  up.  Above  30  Hz  the  steady-state  pulse  energies  drop  rapidly, 
presumably  due  to  heating,  so  that  the  total  laser  power  also  peaks  at  30  Hz.  The  maximum 
average  power  observed  from  this  gas  mixture  was  8  W.  Other  mixtures  gave  peak  average 
powers  of  up  to  18  W  at  50  Hz. 

Figure  16  shows  the  effects  of  varying  the  gas  fill  pressure  at  the  constant  repetition  rate 
of  30  Hz.  The  energy  transfer  to  the  discharge  peaks  at  25  Hz  although  the  laser  pulse  energy 
continues  to  increase  up  the  maximum  fill  pressure  of  35  Torr.  Above  35  Torr  the  discharge 
was  erratic  with  high  voltage  breakdown  occurring  at  the  output  terminals  of  the  pulsed  power 
supply. 

Additives  to  the  gas  mixture,  such  as  Xe,  were  tried,  without  much  success,  to  see  if  the 
impedance  of  the  discharge  could  be  significantly  lowered.  Figure  17  shows  the  effect  of  adding 
water  vapor  to  the  gas  mixture  (which,  in  this  case  included  1%  of  Xe).  The  water  vapor  does 
decrease  the  fluctuations  in  the  laser  output  energy  and  decrease  the  slow  rate  of  energy  decay 
observed. 
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Figure  17:  Adding  water  vapor  to  the  discharge  decreases  the  slow  rate 
of  the  decay  in  pulse  energy. 
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APPENDK 


RIk  gridfitzx  Pafo:  1 


•dafln*  NFRINGE  1200 
•d«fln«  NGRID  41 
tdcfin*  NPASS  10 
td*£ln«  NOPTIMIZE  8 
•daflna  NSMOOTH  2 
tdaflna  NPOINTS  22 
tdeflne  ZFRAC  (float) 0.2 
*d«fln«  TOL  (float)5.0«-S 

/*  GrldFltZ.c  VI. 0  22May92  S.F.Fulghum  Scienca  Raaaarch  l.aboratocy 

This  vsrslon  utlllzas  a  tru«  minimization  of  the  surface  rather 
than  an  average. 

This  version  subtracts  the  polynomial  fit  of  piston,  tip  and  tilt 
from  the  INPUT  DATA,  BEFORE  proceeding  with  the  minimization  fit 
and  Initializes  the  grid  with  the  polynomial  fit  without 
piston, tip, tilt. 

Calculates  a  wavefront  from  standard  Mach-Zehnder  fringe  data 
over  an  NGRID  x  NGRID  grid  using  a  minimum  surface  area  constraint 
to  define  the  values  at  grid  points  not  fixed  by  data. 

Fringe  data  is  stored  In  the  structured  array  f[l]  of  type  fringe 
which  contains  the  fringe  x, y, z  data,  the  location  of  the  grid 
In  which  the  fringe  point  la  found  and  the  Interpolation 
weights  which  determine  the  surface  value  In  terms  of  the  values 
at  the  four  corners  of  the  grid,  LL,LR, UL, UR; 

Grid  data  Is  stored  In  the  2-0  array  g[row](col]  of  type  grid 
and  Is  accessed  through  pointers  to  each  row  of  the  array 
In  the  style  that  Numerical  Recipes  uses  for  arrays. 

g[row] [col]  contains  the  x,y  and  z  values  of  the  grid  point, 
a  Boolean  type  as  to  whether  or  not  the  grid  point  Is  Inside  the 
aperture  of  the  data  (l.e.  for  round  or  rectangular  apertures), 
the  number  of  fringe  points  affected  by  changes  In  this  grid  value 
and  a  structured  array  (type  grldfrlnge)  with  Information  about 
the  fringe  points  affected  by  this  grid  point. 

The  type  grldfrlnge  has  two  elements,  the  first  of  which  Is 
the  Integer  Index  of  the  fringe  data  point  In  question  and 
the  second  of  which  Is  the  weight  the  fringe  point  In  question 
plays  In  the  determination  of  the  grid  point  value. 

In  this  program  rows  Increase  from  the  bottom  of  the  aperture 
to  the  top  of  the  aperture  and 

columns  Increase  from  the  left  to  right  so  that 
col«0  ->  X— 1  to  col-NGRID-1  ->  x«+l; 
row-0  ->  y— 1  to  row-NGRID-1  ->  y-+l; 

INPUT  FILES: 

.NFG  Is  the  normalized  fringe  data  output  from  READFNG 
.POL  Is  the  polynomial  fit  coefficients  from  POLTFIT 

OUTPUT  FILES: 

.GPF  Is  the  polynomial  coefficients  (without  piston,  tip,  and  tilt) 
as  calculated  at  the  grid  positions. 

.GPT  Is  the  polynomial  coefficients  Including  tip  and  tilt 
evaluated  at  the  grid  positions. 

.GMF  Is  the  minimization  fit  to  the  original  data  (from  which 
the  original  tip  and  tilt  had  been  subtracted) .  Note  that 
any  additional  tip  and  tilt  resulting  from  the  minimization 
remains  In  this  data,  (possibly  should  be  changed) 

.GMT  Is  the  minimization  fit  with  the  original  tip  and  tilt 
added  bach  in  so  that  fringes  can  be  calculated. 

.GHD  Is  a  difference  flit  between  the  .GMT  file  and 
a  reference  .GMT  file. 


The  .GMT  output  file  contains  the  minimization  fit  (to  data  from  which 
piston,  tip  and  tilt  had  )3een  subtracted)  with  the  original  tip  and 
tilt  put  back  In.  If  GRIDFITS  Is  run  on  a  null  shot.  It  will  generate 
a  .GMT  file  suitable  for  use  as  a  reference  surface. 


The  .GMD  file  la  the  DIFFERENCE  between  a  .GMT  file  and  a  REFERENCE  .GMT 
file  so  that  an  added  aberration  can  be  dstermlned  which  INCLUDES 
any  added  tip  and  tilt. 
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*/ 

tdafln*  TRUE  1 
•define  FALSE  0 
•include  <atdlo.h> 

•Include  <pcocese.h> 

•  include  <inelloc.h> 

•include  <meth.h> 

•include  <strlng.h> 

/*  •**•*•«*•****•**•****••*>****•••••••**•*•**•«**«•*••••*••*•***••*•**  «/ 

/*  function  declarations  */ 

float  polydnt  m,  float  x.  float  y)  i 

float  polyflt (float  *pcoef,  float  x,  float  y) ; 

float  RHSflt (void) ; 

void  movo_gr id  point (void) ; 

void  smooth_9rId(vold) ; 

void  opt  iml ze_gr id  ( vo id )  ; 

float  L(float~zl; 

void  zllmlts (float  *zmln,  float  *zmax,  float  *zavg) ; 
void  minimize  surface_tenslon (void) ; 

float  sfbrentTfloat  ax,  float  bx,  float  cx,  float  (*f) (float  x) , 
float  tol,  float  *xmin,  Int  *niter) ; 


/*  Global  Variables 

These  structures  contain  the  input  fringe  positions  and  orders, 
the  grid  values,  the  Information  as  to  which  fringe  points  are 
affected  by  each  grid  value  and  the  constants  required  for 
interpolating  the  surface  between  the  grid  points  at  the  positions 
of  each  of  the  input  fringe  data  points. 


•/ 


struct  fringe 

( 

float  x;  /•  X  value  of  fringe  point  •/ 

float  y;  /•  y  value  of  fringe  point  •/ 

float  z;  /*  order  of  fringe  (Integral  but  float  for  calculation)  */ 
int  rLL;  /*  row  of  LL  corner  of  grid  square  containing  the  point  */ 

int  cLL;  /•  column  of  lower  left  corner  •/ 

float  wLL;  /*  weight  of  LL  grid  value  for  the  x, y  value  of  the  point  */ 
float  wLR; 
float  wUL; 
float  wUR; 


struct  grid 


( 


float  xg; 
float  yg; 
float  zg; 
int  in; 
int  nf; 


/*  x  value  at  the  grid  point 
/*  y  value  at  the  grid  point 
/*  z  value  at  the  grid  point 
/*  1  if  Inside  the  aperture. 


*/ 

*/ 

(the  final  output) 
0  if  outside  */ 


/*  number  of  fringe  points  affecting  grid  point 


int  fi(  NPOINTS  ];  /*  vector  of  fringe  Indices  */ 


*/ 

*/ 


struct  fringe  *f; 
int  npts; 
struct  grid  **g; 
float  **temp; 
int  row, col; 

void  malnO 

{ 

FILE  *fp; 

char  basename [40] ,name [40] , inf lie [80] , polyfile [80] ,polyf itf lie [80] ; 

char  polytipflle [80] ,  mlnfltflle [80] , mintlpf lie [80] ; 

char  subflle[80] ,dlfflle[80] ; 

int  i,  j,rLL,cLL,nf,  lfp,maxnf; 

int  ipass, lopt, ismooth; 

slze_t  vecslze; 

float  delta, xg, yg, xmin, ymin,  xmax, ymax, f rLL, fcLL, u, v, wfp; 

float  zero- (float) 0; 

float  two-(float)2.0; 

float  one-(float) 1.0; 

float  sum, dlff ,pcoef [16] ; 
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floac  nu, rmspol, rmsO; 
float  avgdlffO.avgdlff; 
float  xIX.  cut,  zUL,  «UR; 
float  wLL, wLK, hUL, wUR; 
float  gsub; 


/. - - - 

/•  get  filenames  •/ 

prlntf ("enter  normalized  fringe  filename,  d:\\frlnge\\dat\\*.IIFG  assumad\n*) 
strcpy  (basename,  "drWfrlngeWdatW") ; 
gets (name) ; 

strcat (basename, name) ; 
printf ("\n") ; 

/*  Input  data  file  of  normalized  fringes  */ 
strcpy (Inflle, basename) ; 
strcat (Inflle, " .NFG") ; 
printf ("Inflle:  %s\n", Inf lie) ; 

/*  Input  data  file  of  polynomial  fit  coefficients  "/ 
strcpy (polyf lie, basename) ; 
strcat (polyflle, " .POL") ; 
printf ("polyflle:  %s\n", polyf lie) ; 

/•  output  data  file  of  Initial  polynomial  fit,  •/ 

/*  WITHOUT  tile  piston,  tip  and  tilt  coefficients  •/ 
strcpy (polyf Itf lie, basename) ; 
strcat (polyf Itflle, ".GPF") ; 
printf ("polyfltflle:  %a\n",polyfltfllo) ; 

/*  output  data  file  of  Initial  polynomial  fit,  WITH  original  tip  and  tilt  •/ 
strcpy (polytlpf lie, basename) ; 
strcat (polytlpf lie, " .GPT") ; 
printf ("polytlpflle:  %s\n",polytlpflle) ; 

/•  output  data  file  of  minimization  fit,  WITHOUT  original  tip  and  tilt  •/ 
strcpy (mlnfltflle, basename) ; 
strcat (mlnfltflle, ■ .GMF") ; 
printf ("mlnfltflle:  %s\n", mlnfltflle) } 

/*  output  data  file  of  minimization  fit,  WITH  original  tip  and  tilt  */ 

/*  for  plotting  contours  to  chec)i  against  Input  fringes  */ 
strcpy (mlntlpf lie, basename) ; 

St  Tcat (mlntlpf lie, ".GMT") ; 

printf ("mlntlpf lie :  %s\n", mlntlpf lie) ; 

/*  output  data  file  of  difference  (between  two  minimization  fits  */ 

/•  without  tip  and  tilt  on  either  •/ 
strcpy (dlf file, basename) ; 
strcat (dlfflle, " .GMD") ; 
printf ("dlfflle:  %s\n", dlfflle) ; 

/*  Input  data  file  of  minimization  fit  to  SUBTRACT  from  .GMT*/ 

/*  Including  tip  and  tilt  */ 

printf ("enter  subtraction  filename,  .GMT  assumedNn") ; 

/*  gets (basename) ;  */ 

strcpy (basename, "d: \\frlnge\\dat\\cln0S024") ; 

strcpy (subfile, basename) ; 

strcat (subfile, " .GMT") ; 

printf ("subfile:  %s\n\n", subfile) ; 

printf  ("hit  Icey  to  continue\n")  ; 
getchar  0 ; 


/♦ - 

/*  allocate  memory 

NFRINGB  is  the  number  of  fringe  points  that  can  be  stored  (typ.  1000)  . 
NGRID  is  the  size  of  the  grid,  an  odd  value  to  insure  a  grid  point 
at  the  exact  center  of  the  grid  (typ.  41) . 

NPOINTS  Is  the  maximum  number  of  Input  fringe  points  that  can  be 
assigned  to  a  given  grid  point.  The  value  required  will  vary  with 
the  density  of  sampling  the  input  fringes  and  the  density  of  grid 
points  (typ.  10) . 


*/ 


vecslze- (slze_t)  NFRINGE  "slzeof (struct  fringe) 
f> (struct  fringe  *) malloc (vecslze) ; 
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t  (f->NULL)  ( 

prlncf ("error  elloceclng  apece  for  fringe  data\n"); 
exit  (1)  ; 


vecslze**  (slze^t)  NGRID  "slzeof  (atruct  grid  *)  ; 
g-(atrucC  grla  **)ffialloc(vecalze) ; 

If  (g— NULL)  { 

prlncf ("error  allocating  apace  for  grid  row  polnteraVn") ; 
exit(l); 

) 

vecalze- (alze_t)  NGRID  "alzeof (atruct  grid); 
for(l-0;l<  NGRID  ;1++)  1 

g[l]-(struct  grid  *) malloc (vecalze) ; 

If  (gdl— NOLL)  ( 

prlntf ("error  allocating  apace  for  grid  row  %d  data\n",l); 
exlt(l); 

> 


vecalzea (alze  t)  NGRID  "alzeof (float  •) ; 
temp-(float  *'^)  malloc  (vecalze) ; 

If  (temp— NOLL)  ( 

prlntf ("error  allocating  space  for  temp  grld\n") ; 
exit  (1)  ; 

) 

vecalze- (alze_t)  NGRID  "alzeof (float) ; 
for(l-0;l<  NGRID  ;1++)  { 

temp [1]- (float  *) malloc (vecalze) ; 

If  ( temp [1]— NULL)  { 

prlntf ("error  allocating  apace  for  temp  grid  row  %d\n"<l); 
exit (1) ; 

) 


/* - ./ 

/"  Initialize  fringe  and  grid  arrays  "/ 

for (1-0; i<  NFRINGE  ;1++)  ( 
f [1] .x-zero; 
fill .y-zero; 
f (1] .z-zero; 
f [1] .rLL-0; 
f (11 .cLL-0; 
f [ij .wLL-zero; 
f  [1]  .wLR  'Zero; 
f [1] .wUL-zero; 
f (ij .wUR-zero; 


delta  -  two/ ( (float) NGRID  -one); 
yg- (float) -1; 

for (row-0 ;row<  NGRID  ;row++)  ( 
xg- (float) -1; 

for (col-0;col<  NGRID  ;col++)  { 
g(row] (col] ,xg-xg; 
girow] (col] .yg-yg; 
g(rowj (col] .zg-zero; 
girow] [col] .in-  TRlffi  ; 
girow] icol] .nf-0; 

for(i-0;i<  NPOINTS  ;i++)  g(row] (col] .fi (i]-0; 
xg+-dolta; 

) 

yg+-delta; 


/*  -  */ 

/"  read  polynomial  coefficients  */ 

prlntf ("reading  file  of  polynomial  coefflclentsXn") ; 
fp-fopen (polyfile,  "r") ; 

If  (fp— NULL)  { 

prlntf ("error  opening  polynomial  file:  %s\n",polyflle) ; 
exit(l); 

) 

for(i-l;i<-15;i++)  { 

fscanf(fp, "  %f", s (pcoef (1] ) ) ; 
prlntf i"%d  %g\n", 1, pcoef (1] ) ; 

} 

fclose (fp) ; 
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/• -  ./ 

/*  c«ad  frlng*  (l«t« 

This  ssctlon  rssds  sn  arbitrarily  long  list  at  data  trlplsts: 
x-posltlon  y-posltlon  frlng*  order 


prlntf ("reading  fringe  data\n“) ; 
fp-fopen (Infile, "r") ; 

If  (fp— NOLL)  ( 

prlntf (“error  opening  fringe  data  file;  %s\n“, Inflle) ; 
exit(l) ; 

) 

i“0; 

while (! feof (fp) )  { 

f3canf(fp.“  %f  %f  %f“,t(f[l].x),4(f[ll.y).t(f(l].z)); 
i++; 

) 

npts-1-1; 

If  (ferror(fp))  ( 

prlntf (“error  reading  fringe  data  flle\n“); 

fclose (fp) ; 

exlt(l); 

) 

fclose (fp) ; 

prlntf (“%d  fringe  points  read  from  file\n“,npts) ; 


/* -  ./ 

/*  subtract  piston,  tip,  tilt  from  fringe  data 

pcoef[l]  Is  constant,  pcaef[2]  Is  x,  pcoef[3]  Is  y 

•/ 

prlntf (“subtracting  piston,  tip,  tilt  from  fringe  data\n“); 
for (i-0;l<npts;i++)  { 
sum- (float) 0; 
for(j-l; j<-3; j++)  { 

sum+-pcoef [ j) “poly ( j, f [1] .x, f [1] .y) ; 

) 

f  [1]  .z— sum; 

) 


/*  -  ./ 

/•  Initialize  grid  with  polynomial  fit 
WITHOUT  piston,  tip  and  tilt. 

•/ 

for (row-0, 'row<  NGRID  ;row++)  ( 
for (col*-');col<  NGRID  ;col++)  { 
sum-zero; 

for (1-4; i<-15; i++)  ( 

sum+-pcoef [1] “poly (1, 9 [row] [col] .xg,g [row] [col] .yg) ; 

) 

g[row] [col] .zg-sum; 

) 

) 


/. -  ./ 

/*  write  Initial  polynomial  grid  data  to  .GRT  file 

and  the  fit  WITHOUT  piston,  tip  and  tilt  to  .GPF  file 
Uses  exponential  notation  to  comply  with  Mathematlca 
for  plotting  results. 


fp-fopen (polyf itf lie,  “w“)  ; 
if  (fp— NULL)  ( 

prlntf (“error  opening  grid  polynomial  file:  %s\n“,polyfltflle) ; 
exit  (1)  ; 

) 

for (row-0 ;row<  NGRID  ;row  ++)  ( 
for (col-0;col<  NGRID  ;col++)  ( 
sum-zero; 

for ( 1-4 ; i<-15; i++) 

sum+-pcoef [1] “poly (l,g [row] [col] .xg,g[row] [col] .yg) ; 
fprlntf (fp, “*-11.3o\n“,sum) ; 

) 

} 

if  (ferror(fp))  ( 

prlntf (“error  writing  grid  polynomial  fit  fHe\n“); 
fclose (fp) ; 
exit (1) ; 

) 

fclose (fp) ; 
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fp-fop«n (polyclpf lie,  "w")  ; 
if  (fp— NULL)  { 

printf (“error  opening  grid  polynomial  file;  %s\n“,polytipfile) ; 
exit (1)  ; 

) 

for (row-0; row<  NGRID  ;row  ++)  { 
for (col-0;col<  NGRIO  ;col++)  ( 
sum-rero; 

for (i-1; i<-15; i++) 

sum+-pcoef [i] »poly (i,g[row] (col] .xg.g[row] (col] .yg) ; 
fprlntf  (fp,  *1-11 .3e\n'',  sum)  ; 

) 

] 

if  (farror(fp))  ( 

prlntf (’error  wricing  grid  polynomial  fic  file\n*); 
fclose (fp) ; 
exit  (1) ; 

1 

fclose (fp) ; 


/.  - - -  ./ 

/*  distribuCe  fringe  data  points  among  the  grid  points 

This  section  looics  at  each  of  the  input  fringe  data  points 
and  determines  which  four  grid  points  surround  it. 

It  then  writes  this  Information  into  the  data  structures  for 
each  of  these  four  grid  points  along  with  the  interpolation 
constants  (weights)  for  the  fringe  point. 

The  interpolation  scheme  used  for  the  surface  between  four 

grid  points,  lower  left  LL,  lower  right  LR,  upper  left  UL 

and  upper  right  UR  is  as  follows.  The  fractional  distance 

of  a  fringe  data  point  position  from  the  left  edge  cf  the  surface 

towards  the  right  edge  of  the  surface  is  *u*.  The  fractional 

distance  from  the  bottom  edge  to  the  top  edge  of  the  surface  is  *v*. 

The  suface  value  *z”  is  defined  by 

z-  z_LL(l-u)  (1-v)  +z_UL(l-u)  (v)  4'Z_LR  (u)  (1-v)  +z_UR(u)  (v)  . 

This  surface  is  continuous  for  all~values  of  the  grid  points 
but  its  derivative  is  not.  For  a  point  cn  the  line  connecting  two 
grid  points  the  interpolation  reduces  to  simply  a  linear  interpolation 
between  only  those  two  points.  Note  that  the  surface  is  not  flat 
so  the  fringes  that  it  predicts  are  curved  within  a  surface. 

The  values  "u"  and  *v"  are  determined  by  calculating  the  float 
value  of  the  grid  position  of  the  fringe  point  (i.e.  row-2.15  col-20.45) 
and  subtracting  the  Integral  portion.  The  integral  portion  is  used  to 
determine  the  index  of  the  lower  left  LL  grid  point  of  the  surface 
in  which  the  fringe  point  is  found.  The  value  (1-u) (1-v)  is,  for 
example,  the  weight  of  the  lower  left  grid  point  in  determining  the 
value  of  the  surface.  It  is  also  the  weight  of  the  fringe  point  when 
determining  the  optimum  value  for  the  lower  left  grid  point  that 
minimizes  the  error  between  the  fringe  predicted  by  the  surface  and  the 
measu’'sd  fringe  position.  If  the  fringe  point  is  very  close  to  the 
upper  right  corner  of  the  surface  (u  approx.  1) (v  approx  1)  then 
it  is  primarily  affected  by  the  upper  right  grid  point 
and  not  the  lower  left  grid  point. 


for  (1-0;  Knpts;  1++)  ( 

fcLL-(f[l]  .x-t-one) /delta;  /•  fringe  position  in  grid  units  •/ 

cLL- (Int) fcLL;  /•  integral  porlon  of  the  position  •/ 

f (i] .cLL-cLL;  /•  write  this  to  fringe  structure  •/ 

frLL- (f (1] .y+one) /delta; 

rLL-(lnt) frLL; 

f(l].rLL-rLL; 

u-fcLL- (float) cLL;  /*  fractional  position  for  weighting  •/ 
v-frLL- (float) rLL; 

f (1] .wLL-(one-u) * (one-v) ;  /*  weighting  values  •/ 

f (i] .wUR-u*v; 
f (ij .wOL-(one-u) *v; 
f [ij .wLR-u* (one-v) ; 

row-rLL;  /*  write  data  to  LL  grid  point  data  structure  */ 

col-cLL; 

wfp-f (1] .wLL; 

if  ({row>-0)  SS  (row<  NGRID  )  tt  (col>-0)  (coK  NGRID  »)  ( 

nf-g(row] [col] .nf;  /*  fringe  point  counter  starts  at  0  */ 

if (nf<  NPOINTS  )  (  /*  index  of  fringe  point  affected  */ 

g[row] (col] .fi [nf]-l; 
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» 

nf+'f;  /*  IncraaanC  fringe  point  counter  */ 

g(row] [col] .nf-nf;  /*  update  fringe  point  count  */ 


row-rLL+1;  /•  OR  grid  point  •/ 

col»cLL+l; 

wfp-f [1] .wOR; 

lf((row>>0)tt(row<  NGRID  )  St  (col>-Q)  tt  (coK  NGRIO  ])  ( 

nf«g[row] [col] .nf;  /•  fringe  point  counter  starts  at  0  •/ 

lf(nf<  NPOINTS  )  (  /•  first  point  stored  in  fi[0]  •/ 

g[row] [col] .fi [nf]"i;  /*  Index  of  fringe  point  •/ 

} 

nf4"f;  /*  Increment  fringe  point  counter  */ 

g[row] [col] .nf-nf;  /*  update  fringe  point  count  */ 


row-rLL+1;  /*  OL  grid  point  •/ 

col-cLL; 

wfp-f [1] .wUL; 

If  ((row>-0)st(row<  NGRID  )  tt  (col>-0)  ts  (coK  NGRIO  )) 


nf-g[row] [col] .nf; 
if(nf<  NPOINTS  )  ( 
g[row] [col] .fl[nf]-i; 

} 

nf++; 

g[row] [col] .nf-nf; 

) 


[ 


/*  fringe  point  counter  starts  at  0  */ 
/*  point  1  stored  In  fl[0]  •/ 

Index  of  fringe  point  */ 


Increment  fringe  point  counter 
update  fringe  point  count  */ 


row-rLL;  /•  LR  grid  point  •/ 

col-cLL+l; 

wfp-f [1] .wLR; 

lf((row>-0)SS(row<  NGRID  )  SS  (col>-0)  SS  (coK  NGRID  )) 


nf-g[row] [col] .nf; 
if(nf<  NPOINTS  )  ( 
g[row] [col] .fi[nf]-i; 

) 

nf++; 

g[row] [col] .nf-nf; 

) 


( 


/*  fringe  point  counter  starts  at  0 
/*  point  1  stored  in  fi[0]  •/ 

/»  Index  of  fringe  point  •/ 


Increment  fringe  point  c  >anter 
update  fringe  point  count  •/ 


'/ 


)  /•  end  of  fringe  point  loop  •/ 

printf("%d  fringe  points  dlstrlbuted\n",npts) ; 

maxnf-0; 

for (row-0;row<NGRID;row++)  ( 
for (col-0; col<NGRIO;col++)  ( 

if (g(row] [col] .nf>maxnf)  maxnf-g[row] [col] .nf; 

) 

) 

printf ("maxnf-%d,  space  allocated  for  %d\n",maxnf,  (int)NPOINTS  ); 
if(maxnf>-  NPOINTS  )  ( 

prlntfCtoo  little  space  allocated  for  frlngesXn”) ; 
exit(l); 

} 


/. - ./ 

/*  calculate  Initial  fit  error  */ 

rmsO-RMSflt 0 ; 

prlntf ("Initial  RMS  fit  error-  %g\n",rmsO); 

/*  «**«tt««»»*«««tft*«*«**«t**«**»t«*«****#f»»«««««**«f««*t*****if**««t*  */ 

for  (lpass-l;ipass<-  NPASS  ;lpass+-t')  ( 
prlntf ("pass  %d\n", Ipass) ; 

for (lsmooth-l;lsfflooth<-  NSMOOTH  ;lsmooth+>)  smooth_grld () ; 
rms-RMSflt () ;  ~ 

prlntf ("  after  smooth  rms  error-  %g\n",rm3); 

for (lopt-l;iopt<«  NOPTIMIZE  ;iopt++)  optlmlze_grid() ; 
rms-RMSflt  0; 

prlntf ("  after  opt  rms  error-  %g\n\n",rms) ; 

)  /•  end  optimization  loops  */ 

/*  «f«t»*«***«**«t««»«***t*«*««f**fl**i«f*lf«***f**t«l*««*«tl«*t***ttt*t  */ 

/* - */ 

/*  calculate  final  fit  error  •/ 
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prlntf("\n  Initial  RMS  fit  error-  %g\n",nnaO); 
prlntf(*\n  final  RMS  fit  error-  agNo'^nna); 


/. - 

/*  write  final  grid  data  to  .GMF  file 

Uses  exponential  notation  to  comply  with  Hathematlca 
for  plotting  results. 


prlntf ("writing  grid  minimization  fit  In  .GMF  flle\n") ; 
fp-fopen (mlnf Itflle, "w") ; 

If  (fp— NULL)  ( 

prlntf ("error  opening  minimization  grid  fit  file:  %s\n", mlnf Itflle) ; 
exit (1) ; 

) 

for (row-0 ;row<  NGRID  ;row  ++)  ( 
for (col-0 ;col<  NGRIO  ;col++)  ( 

fprlntf (fp,"%-11.3e\n",g[row] (col] .zg); 

} 

} 

If  (fertor(fp))  ( 

prlntf ("error  writing  grid  minimization  fit  flle\n") ; 
fclose (fp) ; 
exit (1) ; 

) 

fclose (fp) ; 


- - - */ 

/*  write  final  grid  data  to  .GMT  file 

In  which  the  polynomial  piston,  tip,  tilt  are  put  bac)c  In 


prlntf ("writing  minimization  fit  with  original  tip  and  tilt  In  .GMT  flle\n") 
fp-fopen (mlntlpflle,  "w") ; 

If  (fp— NULL)  { 

prlntf ("error  opening  fit  minus  polynomial  file:  %s\n", mlntlpflle) ; 
exit(l); 

} 

for (row-0 ;row<  NGRID  ;row  ++)  ( 
for {col-0;col<  NGRID  ;col++)  { 
sum-zero; 
for (i-l;l<-3;l++) 

sum+-pcoef (l]*poly(l,g(row] (col] .xg,g[row] [col] .yg) ; 
fprlntf (fp, ■%-11.3e\n", (g(row] (col] .zg+sum) ) ; 

) 

If  (ferror(fp))  { 

prlntf ("error  writing  fit  minus  polynomial  file:  %s\n", mlntlpflle) ; 

fclose (fp) ; 

exit(l); 

) 

fclose (fp) ; 


/* - */ 

/*  write  final  grid  data  to  .GHD  file 

Uses  exponential  notation  to  comply  with  Mathematlca. 

Subtract  a  set  of  grid  points  first  using  a  full  file  name. 


prlntf ("reading  subtraction  file:  %s\n", subfile) ; 
fp-fopen (subfile,  "r") ; 

If  (fp— NULL)  ( 

prlntf ("error  opening  subtraction  file:  %s\n”, subfile) ; 
exit(l); 

] 

for (row-0 ;row<  NGRID  ;row  ++)  { 
for (col-0;col<  NGRID  ;col'f'«')  ( 

fscanf (fp, "  %e",C (temp (row] [col])); 

) 

} 

If  (ferror(fp))  ( 

prlntf ("error  reading  subtraction  file:  %s\n*, subfile) ; 
fclose (fp) ; 
exit  (1)  ; 

) 

fclose (fp) ; 

prlntf ("writing  difference  file:  %s\n",dlfflle) ; 
fp-fopen (dlfflle, "w") ; 
if  (fp— NULL)  { 
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printf ("error  opening  grid  difference  file:  %a\n",dlfflle) ; 
exlt(l); 

) 

for (row«0;row<  NGRIO  ;ro»  ++)  ( 
for (col-0;col<  NGRID  ;col'*’>)  ( 
sum-zero; 
for (1-1;1<-3;1++) 

pcoef  [1]  *poly  {l,g[row]  [col]  .xg,g[row]  [col]  .yg)  ; 
fprlntf  (fp.  "l-ll.ieVn”,  (g[row]  [col]  .zg-fsum  -temp[row]  [col] ) ) ; 

) 

If  (ferror(fp))  ( 

prlntf ("error  writing  grid  difference  file:  %s\n",difflle) ; 
fclose (fp) ; 
exlt(l); 

] 

fclose (fp) ; 

/* - 

/•  free  memory  space  •/ 

for(l-0;l<  NGRID  ;1++)  free(g[l]); 
free(g)  ; 

for(l-0;i<  NGRID  ;i++)  free (temp [1] ) ; 
free (temp) ; 
free (f) ; 

) 

fH  tktlr*******************************************^**********************  * 

float  polydnt  m«  float  x«  float  y) 

( 

switch (m)  ( 

case  1:  return  (float) 1; 
brealc; 

case  2:  return  x; 
brealc; 

case  3:  return  y; 
brealc; 

case  4 :  return  x*x; 
break; 

case  5;  return  x"y; 
break; 

case  6:  return  y*y; 
break; 

case  3;  return  x*x*x; 
break; 

case  8:  return  x*x"y; 
break; 

case  9:  return  x*y*y; 
break; 

case  10:  return  y*y*y; 
break; 

case  11;  return  x*x*x*x; 
break; 

case  12:  return  x*x*x*y; 
break; 

case  13:  return  x*x*y*y; 
break; 

case  14;  return  x*y*y"y; 
break; 

case  15:  return  y*y*y*y; 
break; 

default:  prlntf("f  subscript  out  of  bounds \n" ) ; 
break; 

] 


float  polyflt (float  "pcoef,  float  x,  float  y) 

( 

Int  m; 
float  sum; 

sum- (float) 0; 

for (m-l;m<-15;ra++)  ( 

sum-t— pcoef  [m]  "poly  (m,  x,  y)  ; 

) 

return (sum) ; 


float  RMSf it (void) 

[ 


Int  l,rLL,cLL, 
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double  sum, zdlff, zvel; 

sum- (double) 0; 
for (l-0;l<npts;l++)  { 
rLL-f (i] .rLL; 
cLL-f [1] .cLL; 

zvsl-  (double)g[rLL  ] [cLL  ]. zg* (double) £ [1] .wLL 
+ (double) g [rLL+1] [cLL  ] .zg* (double) f [1] .wUL 
-f  (double) girLL+lj  [cLL+lj  .zg*  (double)  £(11  .wOR 
> (double) g( rLL  ] [cLL+l] .zg* (double) £[1] .wLR; 
zdl££- (double) £[1] .z-zval; 
sum'f<>zdlf£*zdl£f  ; 

) 

return  (float) sqrt(  sum/ (double) npts  ); 


void  move_grid_point (void) 

Int  1, l£p, rLL, cLL; 

float  avgdlff; 

avgdlff- (float) 0; 

for (i-0;l<g(row] [col] .nf;i++)  { 

ifp-g[row] (col] .fi [1] ;  /•  Index  of  fringe  point  */ 

rLL-f [ifp] .rLL; 
cLL-f [ifp] .cLL; 
avgdiff+-  f[ifpl.z 

-g[rLL  ] (cLL  ] .zg*f [ifp] .wLL 
-girLL  ] [cLL+1] .zg*f (ifp] .wLR 
-g(rLL+l] (cLL  ] .zg*f (ifp] .wOL 
-g[rLL+l] [cLL+1] .zg*f (ifp] .wUR; 

] 

avgdiff/- (float)g(row] (col] .nf; 

g(row] [col] .zg+-  ZFRAC  *avgdiff; 


/*  smooth  grid  */ 

void  smooth  grid (void) 

{ 

float  zsum, nsum; 
prlntf("smo  *); 

for  (row-0;  row<  NGRID;  row+'f)  ( 
for (col-0;  col<  NGRID  ;col++)  ( 

zsum-g[row] [col] .zg+g(row] (col] .zg; 
nsum- (float) 2; 
if(row>0)  ( 

zsum+-g[ row-1] [col] .zg; 
nsum-t—  (float)  1; 

} 

lf(row<  NGRID  -1)  ( 

zsum+-g(row+l] [col] .zg; 
nsum-t— (float)  1; 

) 

if(col>0)  { 

zsum-f-g(row]  (col-1]  .zg; 
nsum-t— (float)  1; 

} 

if(ool<  NGRID  -1)  { 

zsum-i— g[row]  [col-i-l]  .zg; 
nsum-t—  (float)  1; 

} 

temp [row] [col]-zsum/nsum; 

) 

] 

for  (row-0 ;row<  NGRID;  row-f-t-) 
for  (col-0;  coK  NGRID  ;col-t-t-) 
g[row] [col] .zg-temp[row] [col]; 


/*  optimize  grid  •/ 

void  optimize  grid (void) 

( 

printf("opt  "); 

for  (row-0 ;row<  NGRID;  row-f-l-)  { 
for  (col-0;  col<  NGRID  ;col-t-t-)  { 

if  (g[row]  [col]  .nf— 0)  mlnimlze_sur£ace_tenslon()  ; 
else  move_grid_point 0 ;  ~ 
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) 

) 

) 

/*  Langth  of  Grid  Connactlona,  L  (actually  raturnad  as  L^2)  */ 

float  L(float  z) 

( 

float  zdlf.Ltotal; 

Ltotal-*  (float)  0; 
lf(cow>0)  { 

zdif-z-girow-l] [col] .zg; 

Ltotal-i->sqrt  (  (float)  1  4'zdif*zdlf  ); 

) 

if(row<  NGRID  -1)  { 

zdlf-z-g(row+l] [col]  .zg; 

Ltotal+“sqrt (  (float) 1  >zdlf*zdlf  ); 

} 

lf(col<  NGRID  -1)  { 

zdif»2-g[row] [col+1] .zg; 

Ltotal+“sqrt (  (float) 1  +zdlf*2dif  ); 

) 

lf(col>0)  ( 

zdif-z-g(row] [col-1] .zg; 

Ltotal>-sqrt  (  (float)  1  'fzdlf*zdlf  ); 

) 

raturn  Ltotal; 

) 

/ •  ./ 

/*  Minimum  and  Maximum  of  surrounding  points  */ 

void  zllmlts (float  ■zmln,  float  *zmax,  float  *zavg) 

{ 

float  z,zsum<naum; 


If  (row»>20)  ( 

printf (•row-20\n") ; 

) 

*/ 

*zmln- (float) l.OelO; 

*zmax> (float) -l.OelO; 
zsum> (float) 0; 
nsuma (float) 0; 
if(row>0)  { 

z«g[row-l] [col] .zg; 
zsum+«z; 
nsum-Kx  (float)  1; 
if  (*zmax<z)  *zmax»z; 
if  (*zmin>z)  *zmin-z; 

) 

lf(row<  NGRID  -1)  ( 
z«g[row+l] [col] .zg; 
zsum+xz; 
nsum>- (float) 1; 

If  (*zmax<z)  *zmax-z; 
if  (•zmin>z)  ‘zmin-z; 

) 

if(col>0)  ( 

z-g[row] [col-1] .zg; 

zsum+-z; 

n3um+“ (float) 1; 

If  (*zmax<z)  *zmax-z; 
if  (*zmin>z)  *zmin-z; 

) 

lf(col<  NGRID  -1)  { 
z~g[row]  [col-t-1]  .zg; 
zsum-t->z; 
nsum+- (float) 1; 

If  (*zmax<z)  *zmax-z; 

If  (*zmin>z)  *zmin“z; 

) 

if (nsum> (float) 0)  *zavg-zsum/nsum; 
else  ( 

printf ("nsum-0  logic  errorVn") ; 
exit (1) ; 

} 

} 

/ *  * / 

/*  Minimize  surface  tension  */ 
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void  mlnlmlc*  aurf^c*  tanaion (void) 

( 

Inc  nltar; 

float  z, zfflax, zmln, zavg, znew; 

z-glcow] (col] .zg; 

zllmlca ((zmln,  (zmax, (zavg) ; 

afbcanC (zmin, zavg, zmax, L,  TOL  ,(z,tnlt«r) 

g[cow] [col] .zg-z; 


niK  iMdfiigx  PafK  1 


/*  R««dFnq.c 

r««ds  th«  .FNG  file  containing  th«  fringe  poaltlona  In  plxala 
and  convarta  them  to  a  .NFG  or  Normalized  FrlnGe  file 
conalatlng  of  almple  Hat  of  trlpleta:  x  y  order 
normalized  to  a  -1  0  -fl  aquare. 


•/ 

tlnclude  <atdlo.h> 

•Include  <atrlng.h> 

•Include  <proceaa.h> 

void  malnO 

{ 

FILE  ‘fp, ‘fpln, 'fpout; 

char  baaef llename[80] , name [40] , lnflle(e0] ,outflle[S0] ; 
float  xnorm, ynorm; 

Int  xl,yl,z,xmln,ymln,xwlde,ywlde; 

prlntf (’ReadFNG:  reads  .FNG  file  and  normalizes  poaltlonsXn*) ; 

prlntf(*base  file  name.  d:\\frlnge\\dat\\*.FNG  assumed:"); 

strcpy  (basefllename.  *d:  WfrlngeWdatW")  ; 

scanf ("%s", name) ; 

strcat (basefllename, name) ; 

strcpy (Inflle, basefllename) ; 

strcat (Inflle, " . fng") ; 

strcpy (outflle, basefllename) ; 

strcat (outflle, ".nfg") ; 

prlntf ("reading  file:  %a\n", Inf lie) ; 

prlntf ("writing  file:  %s\n", outflle) ; 

fpln-fopen (Inf lie, "r") j 
lf(fpln— NULL)  ( 

prlntf ("error  opening  Input  fringe  data  file:  %s\n*, Inflle) ; 
exit (1) ; 

) 

fpout“fopen (outflle, "w") ; 

If  (fpout— NULL)  { 

prlntf ("error  opening  output  fringe  data  file:  %s\n", outflle) ; 
exlt(l); 

} 


/» 


*/ 


fscanf (fpln, "  %d  %d  td  td",&xmln,Symln,txwide, tywlde) ; 
fscanf (fpln, "%d  %d  %d", sz,  Sxl, Cyl) ; 
while  (z>-0)  ( 

xnorm-  (float)  -1  (float)  2*  (float)  (xl-xmln)  /  (float)  (xwlde)  ; 
ynorm- (float) 1  - (float) 2* (float) (yl-ymln) / (float) (ywlde) ; 
fptlntf (fpout, "%-6.3f\t%-6.3f\t%d\n", xnorm, ynorm, z)  ; 
fscanf (fpln, "Id  td  Id",  tz,  sxl,  Cyl) ; 

) 

fclose (fpln) ; 
fclose (fpout) ; 

) 

/*  •*•*•*****•*••••••*******••••••*****••*••*•••«•**•****•**••••*•*••**••  • / 
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/•  polyflt.c 

fits  a  wavefront  In  the  form  of  a  power  aeries  In  x, y  to 
a  series  of  fringe  positions 

*/ 

tdeflne  MP  15 
tdeflne  NP  15 
tdeflne  SKIP  0 

tlnclude  <stdlo.h> 
tlnclude  <proceaa.h> 
tlnclude  <malloc.h> 
tlnclude  <math.h> 
tlnclude  <strlng.h> 
tlnclude  "XnrXnrutllh.h" 
tlnclude  "NnrXnr.h" 

/*  procedure  declarations  */ 

float  polydnt  m,  float  x,  float  y) ; 
float  polyflt (float  *p,  float  x,  float  y) ; 

/ .  — ... - .............................................................  «/ 

struct  fringe 

( 

float  x;  /*  X  value  of  fringe  point  */ 
float  y;  /*  y  value  of  fringe  point  •/ 

float  z;  /*  order  of  fringe  (Integral  but  float  for  calculation)  */ 

); 

float  polydnt  m,  float  x,  float  y) ; 

/•  ••***••*•**•••*•««*«*•«**«•••>•**•««••*••■•*••***••••«•••«*•**•****•••  »/ 

void  malnO 

{ 

FILE  ‘fp; 

char  basefllename [40] , name [12] , Inflle [80] , out file [80] ; 

Int  1,  j,  k,  l,m,  n,  npts; 
float  wmax, wmln, "w, ‘b, 'x; 
float  xval,yval, zval; 
float  ••a,**u,»»v; 
float  zero- (float) 0; 
float  sum,dlff, rms; 
struct  fringe  ‘f; 

f" (struct  fringe  *) malloc ( (size  t) 1000*slzeof (struct  fringe)); 

If  (f— NULL)  ( 

prlntf ("error  allocating  fringe  data  spaceXn"); 
exlt(l); 

] 

m-MP; 

n-NP; 

w-vector (1,NP) ; 
x-vector (1,NP) ; 
b-vector (1,MP) ; 
a-matrlx(l/MP,l,NP)  ; 
u-matrlx(l,MP,l,NF) ; 
v-matrlx(l>NP,l,NP) ; 

for(l-0;l<1000;l++)  { 
f [1] .x-zero; 
f [1] .y”zero; 
f [1] .z-zero; 

) 

for ( j-1; J<-MP; j++)  [ 
w[ j]-zeto; 
x( j]-zero; 
b[  j]“zero; 

) 

for(j»l; J<-NP; J++)  ( 
for()c-l;k<-NP;lt++)  ( 
a[j] [k]-zero; 
u(j] [k]»zoro; 
vij] (k]-zero; 

) 

} 
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/. -  ,/ 

/*  read  fringe  data 

Thla  aactlon  reads  an  arbitrarily  long  list  of  data  triplets: 
x-posltlon  y-posltlon  fringe  order 


prlntf ("PolyFlt :  reads  normalized  fringe  positions  .NFG  flleXn*) ; 

prlntf(*  writes  polynomial  coefficient  .POL  flle\n*) ; 

prlntf(”base  file  name,  d: WfrlngeWdatW* .NFG  aasumad:"); 

strepy (basefllename, "d: \\frlnge\\dat\\*) ; 

scanf (*%s", name) ; 

streat (basefllename, name) ; 

strepy (Inf lie, basefllename) ; 

streat (Inf lie, ".nfg“) ; 

strepy (outflle, basefllename) ; 

streat (outflle, " .pol") ; 

prlntf ("reading  file:  %s\n", Inflle) ; 

printf ("writing  file:  %s\n",outflle); 

fp-fopen (Inflle,  "r") ; 
if(fp--NULL)  ( 

printf ("error  opening  fringe  data  file\n"); 
exit (1) ; 

) 

i«0; 

while ( I feof (fp) )  { 

fscanf(fp, "  %f  %f  %f", cxval, tyval, Czval) ; 

f [1] .x-xval; 

fill .y-yval; 

f  [1] .z>zval; 

i++; 

) 

npts“l-l; 
if  (ferror(fp))  { 

printf ("error  reading  fringe  data  flleSn") ; 
fclose (fp) ; 
exit (1) ; 

) 

fclose (fp) ; 

printf ("%d  fringe  points  read  from  file\n",npt8) ; 


/* - ./ 

/*  generate  matrix  and  column  vector  for  minimization, 
b  is  the  column  vector  b_m  «  Sum_l  z_l*poly(m,x  i,y_i). 
a  is  the  matrix  a_ran  ■  Sum_l  poly (m,x_l,y_i) •poTy(n7x_l,y_i) . 

Note  that  a  mn  •  a_nra.  ~  ~  ~  ~  ~ 

Solve  [ [a] ]Tx]-[b]~for  the  coefflcents  [x] . 

•/ 

printf ("generating  matrix  and  vector  for  minlmizatlonNn") ; 
for ( j-1; j<-15; j++)  ( 
sum-zero; 

for (i-0;i<npts;l++)  sum+-f (1] .z»poly ( J, f (ij .x, f (1] .y) ; 
bl jl-sum; 

for  ()t-l;k<-15;k++)  { 
sum-zero; 

for (l-0;i<npt8;l++) 

sum+-poly( j,f [1] .X, f il] .y) "poly (k, f (1] .x,f [1] .y) ; 
a( J] (k]-sum; 

) 

} 


/. - */ 

/*  Solve  [[a]]lx]-[b]  for  the  coefflcents  [x] .  */ 

for(k-l;k<-MP;k++)  for (1-1;1<-NP;1++)  u[k] (l)-alk] (1]; 

printf ("decomposing  matrlxNn") ; 
svdemp  (u,  n,  n,  w,  V)  ; 

/*  find  maximum  singular  value  */ 

wmax- ( f loat ) 0 . 0 ; 

for  0:-l;k<-n;k-t'-f)  If  (w[k]  >  wmax)  wmax-w[k]; 
wrain-wmax* (l.Oe-6) ; 
for  (k-l;k<-n;k++) 

If  (w(k]  <  wmin) { 

printf ("setting  value  %d  of  %g  to  zero\n",k,w[k] ) ; 
w[k]-(float)0.0; 

) 
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prlncf ("begin  beckaubstltutlon\n") ; 
svbkab  (u,  w,  v,  n,  n,  b,  x) ; 

prlntf ("aolutlon  vector\n\n") ; 

for (i-l;i<«NP;l++)  prlntf("%d!  %g\n*. l,x[l] ) > 

prlntf ("Xn") ; 


/•  calculate  fit  error  •/ 


aum>zero; 

for(l-0;l<npta;l++)  { 

dlff-f [1] .z-polyflt (X, f [1] .X, f [1] .y) ; 
auni+-dlff»dlff; 

) 

rma-  (float)  aqrt  ( (double)  aunt)  /  (float)  npta; 
prlntf ("rraa  error-  %g\n",rma); 


/*  write  polynomial  coefflclenta  to  output  file  */ 

fp-fopen (outf lie,  "w") ; 
lf(fp~NULL)  ( 

prlntf ("error  opening  polynomial  data  flleXn"); 
exit  (1) ; 

) 

for  (1-1;  K-NP;  i++)  fprlntf  ( fp,  "%g\n".  x  (1] )  ; 
fcloae (fp) ; 


/*  clean  up  and  go  home  */ 
free (f ) ; 

free_matrlx(a, 1,MP, 1,NP) ; 
freo_matrlx(v, 1,NP, 1,NP) ; 
free_matrlx(u,  1,MP,  1,NP) ; 
freo~vector (b, 1,NP) ; 
free~vector (x, l.NP) ; 
free~vector (w,  1,NP) ; 


float  polydnt  m,  float  x,  float  y) 

( 

switch (m)  ( 

case  1:  return  (float) 1; 
break; 

case  2:  return  x; 
break; 

case  3:  return  y; 
break; 

case  4 :  return  x*x; 
break; 

case  5:  return  x*y; 
break; 

case  6:  return  y*y; 
break; 

case  7 :  return  x*x*x; 
break; 

case  8:  return  x*x*y; 
break; 

case  9:  return  x*y*y; 
break; 

case  10:  return  y*y*y; 
break; 

case  11:  return  x*x*x*x; 
break; 

case  12:  return  x*x*x*y; 
break; 

case  13:  return  x*x*y*y; 
break; 

case  14:  return  x*y*y*y; 
break; 

case  15:  return  y*y*y*y; 
break; 

default:  prlntf ("f  subscript  out  of  boundsXn"); 
break; 

) 
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float  polyflt (float  *p,  float  x 

( 

Int  ms 
float  sun; 

sum> (float) 0; 
for  (m-l;n«»15;ni++)  { 
sum+-p  [ra]  ‘poly  (m,  x,  y)  ; 

} 

return (sum) ; 


,  float  y) 


*  J 


